Interpreting concurrence in terms of covariances in a 
generalized spin star system 
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The quantum dynamics of M pairwise coupled spin i is analyzed and the time evolution of 
the entanglement get established within a prefixed couple of spins is studied. A conceptual and 
quantitative link between the concurrence function and measurable quantities is brought to light 
providing a physical interpretation for the concurrence itself as well as a way to measure it. A 
generalized spin star system is exactly investigated showing that the entanglement accompanying 
its rich dynamics is traceable back to the covariance of appropriate commuting observables of the 
two spins. 
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I. INTRODUCTION 

Quite recently a growing attention has been devoted to 
interacting spin systems jlj, ,2„ ^ also in view of the 
fact that they can be successfully used for gate operations 
in sohd state quantum computation processors IE Sill ■ 
They indeed provide scalable systems that can be eas- 
ily integrated into standard silicon technology. Spin sets 
with an assigned one- or multi-dimensional spatial distri- 
bution appear to be comparatively more promising can- 
didates for the realization of entangled states in matter 
systems. 

Numerous papers published over the last few years wit- 
ness the central role played by spin systems, both from a 
theoretical and applicative point of view, in the emerg- 
ing field of quantum entanglement in solid state physics. 
As an example, and in connection to the subject we are 
going to discuss in the present paper, it is worth explic- 
itly citing those papers dealing with investigation on the 
entanglement get established within a pair of spins be- 
longing to a spin ensemble whose d yna mics is dominated 
by XXX 13, XXZ [13, XYZ [H lH and XY [H [l3 
Heisenberg interaction models. In these cases the analy- 
sis is developed through the evaluation of the concurrence 
function appropriate to estimate the degree of quantum 
correlations get established in a two-spin system. 

It is important to underline that generally speaking 
quantifying and controlling entanglement is a crucial 
challenge both from a theoretical point of view and in 
consideration of its applicative potentialities in various 
field of quantum information. It is in addition important 
to stress that studying spin systems is however of remark- 
able interest in its own. It is well known, for example, 
that a magnetic sample can be suitably analyzed adopt- 
ing a general Hamiltonian model describing a system of 



N spins coupled by exchange interactions with arbitrary 
range and strength 

In this paper we study the quantum dynamics of a sys- 
tem of M Heisenberg pairwise coupled spins | investigat- 
ing in particular the time evolution of the entanglement 
get established within a spin pair subsystem. The anal- 
ysis here reported enable us to bring to light a new con- 
ceptual and quantitative link between the concurrence 
function and an observable quantity having a clear phys- 
ical meaning suggesting, at the same time, a simple way 
of measuring entanglement in a two spin system without 
necessarily reconstructing its state. Our approach fo- 
cusses on a generalized spin star system whose dynamics 
provides an enlightening key to go deep into the physical 
meaning of entanglement and concurrence. In this frame- 
work we succeed in giving a simple recipe to control the 
ability of the system in developing only classical or also 
quantum correlations. 



II. HEISENBERG INTERACTING SPIN 
SYSTEMS 



Our physical system consists of M two-level objects 
whose quantum dynamics is completely described by spin 
operators S^''> = 2ct(') (z = 1, ...,M), a^') being the Pauh 
matrix operator pertaining to the z— th two-level subsys- 
tem. The cartesian components of 5^*' , Sx'' , Sy^^ , and 
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Sz , fulfill the usual angular momentum commutation 
relations and | ± 1^) denote the two eigenstates of a'i^ 
that is ai'^l ± h) = ±1 ± h). 

In this section we wish to keep our presentation gen- 
eral enough to encompass several possible physical sce- 
narios of interest. Thus we adopt the assumption that 
a spin pair belonging to our ensemble of M spins 

{i,j = 1,2,..., M; i ^ j) experiences an Heisenberg 
exchange-like interaction describable as proportional to 
Si^) . The hamiltonian model representing such a 

physical situation in the interaction picture with respect 
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to the free hamiltonian 



M 



(1) 



may be cast in the foUowing general form 

M 

Hi= J2 a^JS^'^ ■ S^^'' 

]>i=l 



(2) 



aij being a real coupling constant left undetermined at 
the moment. Since Hj is a scalar operator with respect 
to S* = S*'*-* and in addition commutes with Hq, then 
the 2— component of S as well as are constants of 
motion. As a consequence, preparing the M spins in an 
eigenstate of Sz, the total system evolves in the corre- 
spondent invariant Hilbert subspace. It is not difficult to 
convince oneself that the reduced density matrix relative 
to a prefixed couple of spins, say the spins i and j, when 
expressed into the standard two-spin basis 



{\h,h 



II,, -1, 



1,; 



-1.)} (3) 



assumes a block diagonal form, each block being biunivo- 
cally singled out by one of the three possible eigenvalues 



of Si"^''= {S^z' + S]?'). Stated another way, the reduced 
density matrix Pij{t), obtained tracing the density op- 
erator of the closed system over the degrees of freedom 
relative to all the M spins except i and j, assumes the 
quite simple form 



.0) 



a{t) 

h{t) c{t) 

c*{t) d{t) 

e(t) 



(4) 



when represented in the ordered basis given by eq. (O 
[l^. We wish to stress that the possibility of writing 
Pij(t) as in eq. Q directly stems from the initial prepa- 
ration being instead independent on the values of the 
coupling constants set {ofy}. On the contrary, the ma- 
trix elements of pij depend on this set too. 

Scope of this section is to prove the existence of dynam- 
ical properties of our system of M spins relying solely on 
the structure of the density matrix pij (t) and not on the 
specific form of the functions a{t), b{t),..., e{t) appearing 
in eq. 

This circumstance appears still more remarkable at the 
light of the fact that, generally speaking, it is not possi- 
ble to exactly solve the dynamics of M interacting spins 
described by the hamiltonian model (01 . Thus results ob- 
tained just exploiting the form of the two-spin reduced 
density operator pij (t) ^ , besides being valid whatever 
the coupling constants are, also provide peculiar tools to 
test approximate dynamical solutions when we are un- 
able to exactly solve the system dynamics. 

With these considerations in mind, we now focus on a 
system of two spins i described at a generic time instant 



f by a density matrix like Q without specifying the an- 
alytic expression of the four population functions a{t), 
b{t),d{t), e{t) and the coherence function c{t). 

Let's begin by observing that in order to guarantee 
that the operator given by equation represents indeed 
a density matrix, the inequality 



(5) 



must be satisfied at any time instant t [l7|. It is easy 
to demonstrate that this relation directly stems from the 
requirement that all the eigenvalues of pij are not nega- 
tive. 

An interesting property assumed by each density ma- 
trix belonging to the class defined by eq. concerns 
the possibility of getting a bridge between the Peres- 
Horodecki (P-H) separability condition 0, ITgl an d the 
measure of entanglement proposed by Wootters j20l | . The 
P-H separability criterium claims that the density matrix 
p of a bipartite system composed by two-level subsys- 
tems, is separable if and only if all the eigenvalues of the 
matrix a obtained from p transposing with respect to the 
indices of only one subsystem, are not negative. In our 
case such a matrix, built up from p^ , may be cast as 
follows 



a{t) c*{t) 

b{t) 

d{t) 

c{t) e{t) 



(6) 



Thus applying the P-H separability criterium, it is easy 
to demonstrate that our two-spin density matrix is 
separable at a fixed time instant t if and only if the con- 
dition 



\cit)\ < y/ait)eit) 



(7) 



is fulfilled. On the other hand, when the density operator 
of two two-level systems has the simple form of cq. (@J) 
it is quite straightforward to evaluate the concurrence 
function Conc{t), introduced by Wootters as a measure 
of entanglement in bipartite system composed by two 
qubits, getting 



Concit) = Max{Q,2{\c{t)\ - ^a{t)e{t))} 



(8) 



Considering eq. © we immediately deduce that at a 
generic time instant t the system is characterized by ab- 
sence of entanglement if and only if condition JT)) is sat- 
isfied in accordance with the P-H separability criterium. 

Let's in addition remark that the presence of entangle- 
ment at a time instant t necessarily implies the existence 
of at least a couple of operators A*^*) and B'^^^ acting on 
the bidimensional Hilbert spaces of the spin i and j re- 
spectively, such that the correlation function 



CAB{t) = Tr{p,,{t)A^^^B^^^} 
-Tn{pMA^'^}Tr,{p,{t)B(^^} 



(9) 
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is different from zero. In eq. © Pk {k = is tfie 
reduced density matrix of the spin k and Trj. denotes 
the trace with respect to its degrees of freedom. At the 
Ught of the results expressed by eq. iQ and ||SJ| a suitable 
pair of operators satisfying eq. © must at least fulfill the 
condition of not being diagonal in the standard basis lOl- 
Thus we are lead to consider the two operators ai and 
Gx''- Exploiting the form of Pij{t) as given by eq. Q it 
is immediate to demonstrate that 

c...At) = mm (10) 

since 

(a«)^rr,{pfc(t)a«} = (fc = ^,J) (11) 

and 

^ Tr{p,,(t)aWa(^-)} = 2^[c(t)] (12) 

If instead we analogously consider the two operators 
Gy^ and Gy^^ we obtain 

a^.Jt) - 23[c(t)] (13) 

A comparison among eqs. (IIOII . 1)1 3|l and ^ brings 
to light the existence of a new direct link between the 
concurrence function and measurable quantities of clear 
physical meaning. We may indeed claim that when differ- 
ent from zero the concurrence function may be expressed 
as 

Condi) = ^fci^Z+ci7. - ^VWW) (14) 

a{t) and e{t) being the probability of finding both the 
spins i and j in their up states or down states respectively. 
It is in addition rather remarkable the fact that if a{t) = 
and/or e(i) = the concurrence as expressed by eq. 
reduces to 

SO that it is zero only if both the two quantum covariances 
Ccr^cr^it) and Cayay{t) vauish. 

This elegant formula suggests in a very transparent 
way that the covariances between the x and y compo- 
nents of the two spins are suitable quantities in order 
to highlight the presence of entanglement in our two- 
spin system. It is however important to emphasize that 
Co-^CT^ ^ [Cay ay 7^ 0) doesn't necessarily imply in its 
own that the system has developed quantum correlations. 
This aspect will appear more clear in the next section 
where we will analyze the exact dynamics of a specific 
system whose hamiltonian model is a particular case of 
that expressed by equation (|2Jl. Concluding this section 
we wish to stress that the quite simple form of Pij{t), 
given by eq. Q), which represents the starting point of 
our analysis, naturally arises in many physical contexts 
not necessarily involving spin i systems. In particular 



the conditions under which eq. (|15|l has been derived 
are verified for example when the dynamics of two iso- 
lated atoms, each in its own Jaynes-Cummings cavity is 
studied pH 123 . |2^ . In this sense, we may claim that 
the concepts and tools reported in this paper are flexi- 
ble enough to be exported into physical situations more 
general than the spin system here envisaged. 

III. GENERALIZED SPIN STAR SYSTEM 

Appropriately choosing the coupling constants a^- the 
quite general hamiltonian model ^ can describe very 
different physical scenarios. Here we wish to focus on 
the case of a system composed by a pair of not directly 
interacting spins each one coupled to every two-level com- 
ponent of a set of spins by a physical mechanism rep- 
resentable by a XY Heisenberg exchange like term. If 
we stipulate the absence of any internal coupling within 
the set of A^ spins, hereafter referred to as the spin bath, 
and indicate by A and B the two spins of the preferred 
couple, also called the central system, the hamiltonian 
model Q assumes the following form 

N N 

HxY = aA{Sx^Sl. + Sy^Sl) + 

N N 
i=l 1=1 

It is worth noticing that specializing the hamiltonian 
model given by eq. ^ into the one expressed by eq. 
H15|l we are supposing that the spin A and B couples 
with any element of the bath with a site-independent 
coupling constant a a or as respectively. 

The interest toward this system, representing a gener- 
alization of those systems reported in literature as spin 
star systems arises from several considerations. 

First of all, as we are going to show, it is exactly solvable 
so that in this case we can explicitly know the functions 
a(t), 6(t),..., e{t) appearing in the density matrix given 
by eq. Q . This circumstance allows us to explore condi- 
tions for the emergence of and to investigate on the time 
evolution of correlations get established within the pair 
A and B. A specific aspect of this study is in particular, 
the possibility of controlling the rising of only classical or 
also quantum correlations between the two not directly 
interacting spins in the central subsystem. Another im- 
portant point worth to be emphasized is the fact that 
under appropriate initial conditions, the concurrence ap- 
pearing and evolving within the A — B subsystem may be 
traced back and interpreted in terms of simple measur- 
able quantities. As a consequence we are in condition to 
suggest a way to measure the concurrence in laboratory 
without being obliged to reconstruct, as usually required, 
the state of the system. 

Let's begin by observing that the possibility of exactly 
solving the dynamics of the system described by Hlt)|l . 
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is strictly related to the existence of other constants of 
motion with respect to the general model 
Introducing the bath collective operators 



N 



j=j:s 



and 



J — Jx ^ iJy 



(17) 



(18) 



and casting the interaction hamiltonian H16|) in the form 

HxY = aA{S+J- + S^J+) + asiS^J- + S^J+) (19) 
it is indeed easy to convince oneself that 

[HxY,J^]^0 (20) 



[HxY, Jrnt] - 



(21) 



J int being an intermediate angular momentum result- 
ing from the coupling of selected at will individual angu- 

— > 

lar momenta of the bath. For example defining J i,j = 
's(i) + ^0), or fc = J*, J + S^C^), we deduce 



[HxY, Jij] — 



[HxY, Jfjh] - 



(22) 



In order to make more evident how the existence of 
all these constants of motion provides the possibility of 
exactly solving the dynamics of the closed system consti- 
tuted by the central one {A, B) and the spin bath, let's 
denote by { | J, M, v)} a. bath coupled basis satisfying 

J^\J,M,v) = J{J +l)\J,M,u) (23) 
J,\J,M,iy) ^ M\J,M,iy) (24) 

u being an integer index taking into account the degen- 
eracy with respect to the quantum numbers J and M. 
As far as the ordered basis of the {A, B) subsystem given 
by eq. ((SJ, we introduce the new more compact notation 



in ,a. 



'in)), 



1,..,4} putting 



k^(l),'^f(l)) 



|U,ls), 
|U,-1b), 
|-lA,li3), 
|-1a,-1b) 



(25) 



We can thus denote by \J,M,i')\a^{n),a^{n))= 
I J, M, V, cr^(n), erf (n)) a generic state of the basis of the 
Hilbert space of the total system obtained as tensorial 
product of the two bases previously introduced. 

It is immediate to convince oneself that the dynamical 
constraints imposed by the constants of motion put into 
evidence at the beginning of this section, allow us to claim 
that, starting from a state |(T^(n), erf (n), J, M, at a 



generic time instant t the state of the system will be at 
most a linear superposition of four states as follows 



\m) = 

it)\af{n),af{n),J,M, i.) + 
Bn.iMit)\ - <jfin), -afin), J, M + <jf{n) + af (n), ly) + 

DnjM{t)\ - af{n),a^{n), J, M + af{n), v) (26) 

We have demonstrated that it is possible to explicitly 
find the exact form of the amplitudes AnjM{t), BnjM{t), 
CnjM{t) and Dnju^t) appearing in eq. (|26|l whatever 
the state of the coupled basis chosen as initial state 
is. Here, for simplicity, we do not give their analyti- 
cal expressions also because, generally speaking, they are 
highly involved. However it is important to stress that, 
knowing the functions A„jm(0: ••■ -DnjA/(0 whatever 
n, J and M are, we, at least in principle, can evaluate the 
temporal evolution of the closed system from a generic 
initial state. 

We now concentrate on a specific initial condition 
namely the one wherein the two central spins are ini- 
tially prepared one up and the other one down whereas 
the bath is in its eigenstate of minimum free energy: 

|V(0)) = |1^,-1^,^,-|,1> (27) 



^ is not 



We remark that the couple J = ^, Af = 
degenerate so that v = 1 only. The time evolution of the 
system starting from this initial condition is character- 
ized by the fact that the probability of finding the two 
spins of interest in the state 1^) is zero at any time 
instant t. This circumstance directly follows from the 
boundaries imposed to the dynamics of the system by 
the conservation of the z component of the total angular 
momentum. Moreover it is possible to prove that in the 
case under scrutiny the matrix element c{t) is real. Thus 
the reduced density matrix of the central system assumes 
the form 



Pi] (t) 





b{t) c{t) 

c{t) d{t) 

e(t) 



where 

m 



(1 



;[cos{2^JN{l + r^^)aAt) + r' 



(28) 



(29) 



z{t) = 



[cos{2^JN{l + r^)aAt) + 



[1 - coB{2^JN{l + r^^)aAt)] 



(30) 



d{t) 



(l + r2)' 



: [1 - cos{2^ N {I + r^)aAt)f (31) 
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FIG. 1: Concurrence function, Cone, against qa^ in correspon- 
dence to N = 100 and r = 1 (point line), r = 10 and r = 10~^ 
(bold line) respectively. 




FIG. 2: Concurrence function. Cone, (bold line), probability of 
finding the state | — 1^, 1^) (point line) and probability of finding 
the state ll'*,— 1^) (continuous line), against aj^t and in corre- 
spondence to N = 100 and r = 1. 



(1 + H) 

In equations H29() - (|32|l r = ^ measures the ratio of the 
two coupUng constants between each component of the 
central system and the bath. Exploiting the results ob- 
tained in the previous section expressed by eqs. H1U|) and 
(|15|l we may thus claim that in the case here analyzed, 
measuring the covariance C^^a^ corresponds to directly 
detect the entanglement arisen between A and B being 

Cone = la.aj = 2|c(i)| (33) 

as easily demonstrable looking at eqs. (|29|l - l|32|l . In par- 
ticular the circumstance that we can explicitly solve the 
dynamics of the system provides the possibility of know- 
ing at a generic time instant t the degree of entanglement 
developed in the central system, starting from a factor- 
ized condition, as measured by the concurrence function 

Figure (1) displays the temporal behaviour of the con- 
currence function Conc(t), obtained putting N = 100, 
in correspondence to three different values of r namely 
r — 10"^, r — 1 and r = 10 respectively. 

As clearly shown, varying the parameter r the con- 
currence function always manifests a periodic oscillatory 
behaviuor immediately deducible taking into account the 
fact that, in the case under scrutiny, Concit) = 2|c(<)|. 
In particular we may state that, whatever r is, there 
exist infinite values of t in correspondence of which 
Concit) = meaning that the central system is separable. 
It is worth emphasizing the remarkable circumstance that 
when 7" = 1 the matrix element c{t) at any time instant 
t assumes its maximum value compatible with relation 
(O, that is \c{t)\ — y/b{t)d{t). In order to verify the 
occurrence of such a saturation of the inequality ij^l at 
any t, it is enough to put r = 1 in eqs. On 
the other hand, it happens that in correspondence to all 
the time instants such that Conc{t) = also the diagonal 
matrix element e{t) vanishes. This means that at these 
time instants all the conditions on the density matrix of 
two spins defining a pure state are verified |l7l | . Thus we 



may claim that when the two spins A and B are equally 
coupled to the spin bath, namely r = 1, the condition 
Concit) — not only reveals absence of entanglement 
in the central system but also guarantees that such two 
spin system is in a pure state. Looking indeed at eqs. 
(|29|l - (|31f) it is immediate to convince oneself that the two 
spins are periodically found in the initial state |1"^, —1"^) 
or in the state | — 1"^, 1^) in which the role played by the 
two spins is exchanged. This behaviour is clearly illus- 
trated in figure (2) where we compare the concurrence 
function with the temporal evolution of the probability 
of finding the state jl'^, —1^) and that of measuring the 
state obtained exchanging the two spins. 

Let's now examine how the system evolves starting 
from another initial condition of experimental interest 
that is 

|V(0)) = 11^,1^,^,-^,1) (34) 

obtained leaving once again the bath in its ground state 
and preparing the preferred couple of spins in the state 
correspondent to the maximum values of Si'^'^\ It is 
possible to prove that under this initial condition and 
assuming that the Hamiltonian model H19|) is invariant 
under the permutation of A with B (a^ = as = a), the 
central system develops correlations when the time goes 
on. More in detail the temporal behaviour of C^^a^ in 
this case is given by 

t^a.a.(i)-^^^sin2c.t (35) 

where 

w = 2V6iV - 4a (36) 

The maximum amount of correlation between and 
is obviously strictly related to the number of external 
spins populating the bath but, however, is larger than 
i as immediately deducible by eq. (|^ . On the other 
hand, evaluating the concurrence function exploiting the 
results presented in section II yields 

a„c(t) = max[0,/w(x)] (37) 
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FIG. 3: Covariance function Ca-^tr^ and concurrence Cone against 
aAt, in correspondence to N = 100 and r = 10 

where 

fN{x) = X + (38) 

2(iV-2)^2iV(iV-l) 

(3iV-2)2 ^ + 

4(1 - N)^/2N{N - 1) + 7V(37V - 2) 
(37V -2)2 

with a; = coswt with w given by eq. (|36|l . It is easy to 
verify that Jn{x) is not positive when a; G [— 1, 1]. Thus 
we may conchide that whatever N is, the concurrence 
function Conc(t) is equal to zero at any time instant t. In 
other words, differently from the case before analyzed, 
under the hypotheses before discussed the concurrence 
function is equal to zero at any time instant t. Stated 
differently the correlations occurring in the central sys- 
tem at a generic time instant t, manifested throughout 
Ccr^a^{t) ^ 0, have only a classical origin and therefore 
we say that such correlations are classical. 

It is interesting to observe that this inability of the sys- 
tem to generate entanglement between the spins A and 
_B when the system is prepared accordingly to eq. H34|) , is 
overcome simply breaking the symmetry condition r — \. 
This fact is clearly shown in figure (3) where we plot the 
two functions Concify and Ca^a^{t) in correspondence to 

= 100 and r = 10. Thus for r ^ 1 the concurrence 
function remains zero for a given interval of time, then 
abruptly increases and once again falls to zero. This be- 



haviour is then periodically recovered reaching maximum 
values of Conc{t) of experimental interest. We wish to 
remark that the evolution of the system in this case is 
characterized by intervals of time during which the sys- 
tem develops only classical correlations {Conc{t) — and 
C„_^„_^ (t) 7^ 0) and intervals of time where quantum cor- 
relations occur {Conc{i) 7^ 0). 



IV. CONCLUSIVE REMARKS 

In this paper we have analyzed the quantum dynamics 
of a system of M Heisenberg coupled spins \ concen- 
trating in particular on a couple of them and looking for 
the time evolution of entanglement developed between 
the two spins. We have demonstrated that it is possible 
to establish a conceptual and quantitative link between 
the concurrence function and easily measurable quanti- 
ties only requiring that at a generic time instant t the 
reduced density matrix describing the two spins of in- 
terest can be put in the simple form given by eq. I@J. 
This circumstance not only allows the possibility of physi- 
cally interpreting the concurrence function introduced by 
Wootters as an useful and powerful mathematical tool to 
estimate the degree of entanglement, but also suggests a 
direct way to measure it. The quantitative relation given 
by eq. 1)14(1 provides indeed the possibility of measuring 
the entanglement without reconstructing the state of the 
two spins. Examining in particular a specific physical sys- 
tem, namely the generalized spin star system discussed 
in section III, we have put into light that detecting the 
covariance function of the two commuting observables 
CT^ and (7^ directly provides the entanglement evolution 
get established between the two central spins. Exploit- 
ing the knowledge of the exact expression of the density 
matrix of the system, pab (t) , starting from an arbitrary 
initial condition, we have disclosed a rich dynamics. In 
particular we have envisaged physical conditions under 
which the two central spins are able to develop classi- 
cal correlations only being in this case Conc(t) = at 
any time instant t. At the same time we have found that 
slightly changing some parameters characterizing the sys- 
tem, quantum correlations appear. 
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